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: + Deseribed and listed herein with concomitant sample input and 
‘output is ‘the - Fortran 1¥-pxogram whitch estimates parameters and standard 
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- @Erorg ‘of estimate per parameters for ee estimated through ~ 
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"1.0 = PROGRAM IDENTIFICATION . es oe 7 3 2 


°2,0 = OBJECTIVE, - -._ -. 


. . \ 
3.0 +,PROGRAM-DESCRIPTION | 


PROGRAM FOR ESTIMATING PARAMETERS .THROUGH MULTIPLE MATRIX SAMPLING WITH 
STANDARD ERRORS OF ESTIMATE APPROXIMATED BY ‘THE JACKKNIFE 


d 


- Program for Estimating Parameters through Multiple Matrix rs ies 
with Standard Errors of Estimate Approximated by the Jackknife/“: — 


5 


“This program estimates selected parameters and standard error of 
estimate associated with each parameter given fultiple matrix 


sampling. “3 


Using item scores sequenced by subtest and, within each subtest, by 
examinee as ‘input, the program computes an estimate of each parameter 
. for each subtest and -then pools results over subtests to provide the 
{ best estimate of each parameter. The parameters estimated within the 
program are: (a) mean test score, (b) second central moment, (c) third 
*_. . central moment, (d).fourth central moment, (e) component_of variance , 
or items, (f)+component of variance for examinees, (g) component of’ 
variance for item x examinee interaction, (h) average item reliability, 
(i) item difficulty for each item tested, and (j) reliability of the ; 
.total score. The equations used to estimate parameters (a) through (d) e 3 
are those given by Lord (1960); the equations for estimating components 
of variance for a completely randomized two-factor design are given: in: 
: maser sources (e.g., Myers, 1966); and the equation for estimating. 
the avérage item reliability is given by Winer (1971). A pooled..{ovéer 
subtests) estimate of each pavameter as well as the estimated standard 
error of estimate associated with.each parameter is accomplished by 
the jackknife procedure, A description of the jackknife procedure and | 
research supporting its use in multiple matrix sampling is given by 
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Shoemaker (1972). The 90, 95, 97.5 and.99 percent confidence intervals 
for parameters (a) through (h) are computed using the t+distribution 
(df = number of subtests minus one) and jackknifed standard errors of 
estimate. The reliability of the total test score is computed by 
projecting the pooled ‘average. item reliability to the total test 

using the Spearman-Brown prophesy formula. Using pooled e.timates of 
moments, the normative relative frequency distribution ts graduated 

by the negative hypergeometric distribution and/or one of the family 


ve 


° 


This program is an improved and expanded version-of an earlier version 
given’ in an appendix of Technical Report No. 34 and abstracted in 

Behavioral Science (Shoemaker, D. M. A Fortran IV program for approxi- , 
mating a normative distribution with the negative hypergeometric 

distribution. Behavioral Science, 197], 16, 414-415.) 
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of Pearson curves ,as opted by the user. The equation for the negative 
hypergeometric distribution is given by. Lord and Novick (1968) and 

the equations for the Pearson curves are given by Elderton and Johnson 
(1969). : 


4.0 - SUBROUTINES AND FUNCTIONS . 


4.1 = Coded Subprograms , “ F 
r 3 . 
, LOAD - Duplication of data base on scratch file number 2; 
No arguments. . 


NEGHGR - Graduation of normative distribution by negative 
hypergeometric distribution. 


. SUBROUTINE NEGHGR (KPOP, XBAR, XVAR) : 


KPOP = Maximum test score 
XBAR ~ Mean test score 
XVAR - Variance of test” Scores 


ITABLE - Calculation of Student's t corresponding to 


prescribed alpha level (two-tailed probability). . 
SUBROUTINE TTABLE (TN ,AREA) 


’ fT - Student's t. ~ , 
N- Degrees, of _ freedom associated with t. 
AREA = Relative area associated with confidence 
interval. 


RDEMT - Subroutine for inputting variable format. 
SUBROUTINE RDFMT (FMT ,N) 


\, - FMT ~ Vector containing format. .. 
N - Number of format cards inputted. 


JKNIFE - Calculation of poolléd estimates ‘of parameters and 
associated standard error of estimate per parameter 
by the jackknife procedure. ‘ 

\ ? 4 . * ° 
SUBROUTINE JKNIFE (U1,U2,U3,U4) _ 
Pooled estimate of mean test score. 

U2 - Pooled estimate.of second central moment. 
U3 ~ Pooled estimate of third central moment. 
Pooled estimate of fourth central moment, 
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SUBTST - Calculation of parameter estimates for each subtest. 
SUBROUTINE SUBTST (NDICH,KPOP,NDISK) | 4 


made ~ Item scoring procedure. 
= dichotomous 
: : = non-dichotomous 
KPOP ~ Maximum test score. 
NDISK = Data base storage. 
O = no storage on scratch file 
1 = data base stored on scratch file 
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" PEARSN - Graduation of normative distribution with one of 
the family of Pearson curves. - ' 


" SUBROUTINE PEARSN (U1,U2,U3,U¢,L1,L2,STEP) 


Ul = Mean test score. ° 
U2 - Second central moment. 
» U3*-e Third central moment. 
U4 - Fourth central moment. 
er Ll - Lower limite of index. 
L2 - Upper limit of index. 
STEP -"Desired step size in going from 
Ll to L2;. : 


POLRT - Computation ‘of the ea and complex roots of a real 
i ceca : . ‘ 


SUBROUTINE POLRT _(XCOF ,M,ROOTR, ROOTI , IER) 


XCOF - Vector of M+ 1 coefficients of the 
. polynomial ordered from “smallest to 
> largest power. 
M ~ Order of polynomial. 


¢ 


ROOTR = Resultant vector of length M Snteiing 
real roots, 

ROOTI - Resultant vector of length M containing 
the corresponding imaginary roots. 


ee IER - Error code return. 
.. 0 = no-error 
‘1 = M less than 1 
2 2 =M greater than 36 
: 3 = unable to determine root with 500 
. iterations on 5 starting values 
4 = High order coefficient is zero ~ oy 
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FRV = Evaluation of cos (x) "e" "dx from -1.570796_to D. 570796 
using Simpson's rule. r 


~ FUNCTION FRV (R,VV) 


= - 6.*(B2-Bl-1.)/(2.*B2-3.%*B1=6.) 


~ SQRT( (R*R* (R-2.)**2*B1)/ (16. *(Rel. ) -BL* (R-2 .)**2)) 


GAMMA - Computation of gamma function using the recursion 
’ ' relation and polynomial SPP vor imetaons 


FUNCTION GAMMA (xX) 
xX - is real argument. . 


4.2 - brary subpeoavene 
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- §.0°= DATA SPECIFICATIONS 
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Sy le Input Fomats 
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Input cena are described in detail at the beginning of 
the program listing. An example of an input data structure 
’ is given in Attachment 1. ‘ 


5.2 = Output. Formats 
Qutput ‘ceantie. from ‘the input. data structure listed in 
sas I is given in Attachment 2. 


6.0 = PROGRAM CONSTRAINTS AND LIMITATIONS ° 


6.1 = Programming Language 


FORTRAN .\ 
6.2 = Vendor | 
a 
University of California at Los Angeles Campus Computing 
> Network. - T 
{ 
1 
v 


6.3 - Storage Requirements 
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cs 
Py 


Total core for execution: . 150K 
See 7.0 for scratch file requirenents. 


* 6.4 = Hardware Configuration 
: ¥ IBM 360 Model 91, Punched cards 


6.5 = Program Parameters 


‘ 


e : t 

6.5.1 - Item scoring: the statistical procedures used within ‘ 
the program assume a uniform scoring procedure for ‘ 
all items in all subtests, ; 


6.5.2 - Parameters estimated: although the program is 
designed primarily to estimate parameters when 
itemg are scored dichotomously (0 = fail, 1 = pass), 
° selected parameters are estimated when items are 
not scored dichotomously--specifically, all parameters 
with the exception of the third and fourth central 
moments. This is a limitation because without 
estimates of these moments it is not possible to 
graduate the normative distribution. Equations for‘ 
the third and fourth central moments when items are | 
: scored non-dichotomously are under development 
: Z currently and will be included in the next revision 
of the program. This limitation in the ‘current 
program is relatively’minor because scoring items 
non-dichotomously is a procedure encountered rarely - 
in PEAGtAees . 


. 


ie , 
6.5.3 = The ee of subtests within any data structure is 
limited to 100. . 


~ 


66 Gu tie nae oe Gens per subtest is limited to 150. 


6.6 - Error Messages 


The Following error messages may be generated by the program: 


© 


i. Errors on‘ card specifying problem set (MAIN) 
(Note: error terminates program) 
2. Error in t-table arg: df= _ | (TTABLE) 
3. Error’in t-table arg: p=_. 7 an (TTABLE). 
BA ag oi ’ : . 7 
4 4. Excessive number of format cards (RDFMT) 
am (Note: error terminates program) 
ee 70 : . ~ ee 


< ‘ ¥ - ¢ A 
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» 5. Invalid limits and/or delta arguments — 
6. Ordinate undefined at this point 

7. M’.LE. 1. for Type XI 

8. Mt out of Fange liga Type VIII 

Me M negative for ‘Type IX 


10, Ganma arg within 0000001 of being a 
‘negative integer 


i Gamma arg 6T 57, eee ans set to 1.E75 


7.0 OPERATING INSTRUCTIONS ; . . 


Submit the program in batch mode in the following fom: 


col 1 * 7 ‘ . 4, va 
* //job card 

// PASSWORD | oO 

// EXEC FORTGCLG, RG=150K 

/ /FORT SYSIN DD # 


main . Pe ; - 


subroutines 
/* 


//GO.FTO2F001 DD DSN=&SCRATCH,UNIT=SYSDA,. 


(PEARSN) 
(PEARSN) . 
(PEARSN) ~ 
(PEARSN)’ 


(PEARSN)* 


(GAMMA) 
(GAMMAY 


// DCB=(RECFM=FBA , BLKSIZE=7280 »LRECL=80) , SPACE= (TRK ,150,RLSE) 


//G0.SYSIN- DD * 
data structures 
/* ; 


// , 


G7 


Input 
parameters 
of data set 


Subroutine \ 
LOAD — 


whe 


/ Subroutine ! 
RDFMT 


Subroutine 
SUBTST 


° 


Subroutine 
- JKNIFE 


Subroutine 
PEARSN ~ 


Subroutine 
NEGHGR 


“9.0 - PROGRAM LISTING 
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x 
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C a . if ‘ t . » 
CE eX ok a kk a ao Ro a kk eo ak A ok ok a tok koi de doe ako ak ak aeak HR KH 
\ 


C \ 
C. 
C ESTIMATIGN OF TEST SCURE PARAMETERS THROUGH MULTIPLE MATRIX SAMPLING 
C WITH STANDARD ERRORS GF ESTIMATE APPROXIMATED BY THE JACKKNIFE 
C PROCEDURE \ 
c : i 
C * \ 
C DAVID M.-SHOEMAKER 1 
C “SOUTHWEST REGIONAL LABORATORY FOR EDUCATIUNAL 
Ci RESEARCH AND DEVELUPMENT : 
€ — 3 
C 
C ORGANIZATION OF DATA STRUCTURES ; - 
: . 
c le DATA CARDS ARE ORGANIZED SEQUENTIALLY INTO PROBLEM SETS. EACH 
C PROBLEM SET CONTAINS (IN THIS ORDER) FHE FOLLOWING CARD SETS: 
c 
Cc, [ TITLE CARD 
C PARAMETER CARD . (' 
DATA BLOCK FOR SUBTEST 1 
DATA BLOCK FOR SUBTEST 2° | . i 
“DATA BLOCK FOR’SUBJEST T° “2 
2. ORGANIZATION GF INDIVIDUAL CARD SETS: . 
‘CARD SET COLUMN "DESCRIPTION 
_ 4 : (INTEGERS RIGHT-JUS TIFIED) 
oe Wg ee a ee ta eee e227 P  8 8 BEE ZF @ Zee ee eer ee Ses ea= Sef = 
Ae TITLE CARD (1 CARD) * 01-72 ALPHANUMERIC TITLE OF 
; PROBLEM SET 
B. PARAMETER CARD (1 CARD) 01-04 - MAXIMUM INTEGER TEST 


SCORE 
05-08 |. NUMBER OF SUBTESTS 
‘ 10... METHOD OF ITEM SCORING 


‘ -* 0=DICHOTOMOUS 
‘ . 1=NON~DICHOTOMOUS 


12 NORMATIVE DISTRIBUTION 


O=NONE DESIRED 
L=NEG.HYPER; 
2=PEARSON CURVE, 
F 3=BOTH 1 AND z 


EE Bara Oe nS NG Om RNC 
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16-30 MEAN . : 
z 31-45 VARIANCE - ° 
. 46-60 3RD CENTRAL MOMENT 
CENTRAL, MOM 


“. 61-75 . 4TH CENTRAL MOMENT ' 
NOTE: THIS PROGRAM MAY ALSO BE 
USED TO GRADUATE NORMATIVE 
DLSTKIBUT IONS WITH THE NECESSARY i 
: MOMENTS CALCULATEO IN ADVANCE « 
4 an IN THIS CASE, THE REQUIRED 
‘ MOMENTS ARE INPUTTED ON COLUMNS 
16 TO 75 (WITH DECIMAL POINT 
PUNCHED GN’ THE CARD). THE NUMBER. ; 
cas OF SUBTESTS (COLUMNS 5 TO 8) te 
SHOULD BE SET EQUAL TO 0. IF 
ESTIMATES OF PARAMETERS ARE TO 
‘BE COMPUTED’ THROUGH MULTIPLE © 
MATRIX SAMPLING, IGNORE COLUMNS 
a 16 THROUGH 75. « 9 
2, a « . 


ANAOAAAKAO 


ee 
4% 


Ce, DATA BLOCK FOR’ SUBTEST I: THE DATA BLOCK FOR EACH SUBTEST 
CONSISTS OF THREE CARD -SETS. 2 


( 


1) PARAMETER CARD (1 CARD) 

2) FORMAT CARD SET (K CARDS) 

3) “DATA CARDS: CONTAINING ITEM SCORES PER EXAMINEE 
- SEQUENCED BY EAA Nee : 


QPOANAANNAGANAGDADANAANAAAAAGIAOA 


THE CRGANIZATION GF THESE CARD SETS IS AS FOLLOWS: 
t / 


PARAMETER CARD 01-04 NUMBER OF EXAMINEES IN SUBGROUP 


Q5-08 NUMBER OF ITEMS PER SUBGROUP 
f ‘ 
: FORMAT CARD: SET 01-72 ' STANDARD FORTRAN IV FORMAT 
_ -'‘ + PUNCHED IN CULUMNS 01-72 ON 
== EACH CARD: AND ENCLOSED IN 
PARENTHESES FOR INPUTTING ITEM 
SCORES FOR EACH EXAMINEE’ IN THE 
: : ITEM-EXAMINEE SAM LE. THE 
NUMBER OF CARDS CONTAINING THE 
FORMAT MAY’ NOT EXCEEU NINE FuR 
EACH ITEM-EXAMINEE DATA SET e 
THE FIRST CARD AFTER THE FURMAT 
CARDS MUST CONTAIN "END UF 
FORMAT" IN CGLUMNS 01-13. 


‘ & 


DATA CARD SbT THE {TEM SCOPES UF EACH EXAMINEE 
; ON SUBTEST SEQUENLED BY EXAMINEL 
- AND ORGANIZED ON CARDS ACCURDING TU 
" FOKMAT SPECIFICATIUN ° 


MOCAAAAIANANAMC ONC. oO 


| ae a 
| ERIC * .-. ; a 
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3e THE PROGRAM WILL PROCESS REPEATED, PROBLEM SETS. 


AOOAYN 
, 
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4, RESTRICTIONS: f8 . oe - 
A. “NUMBER OF SUBTESTS LIMITED TO 100 an 
B. NUMBER OF ITEMS PER SUBTEST-LIMITED TO 150 
Ce WHEN ITEMS ARE SCORED NON-DICHOTOMOUSLY, ONLY THE MEAN TEST 
‘SCORE, VARIANCE OF. TEST SCORES, VARIANCE OF THE ITEM 
. DIFFICULTY INDICES AND TEST RELIABILITY ARE ESTIMATED. 
0. WHEN ITEMS ARE SCORED GICHOTOMOUSLY, THE MEAN TEST SCORE): 
- VARIANCE OF TEST SCORESs 3RD AND 4TH CENTRAL TEST .SCORE 
MOMENTS, VARIANCE OF ITEM OLFFICULTY INDICES AND TEST 
RELIABILITY ARE ESTIMATED. . 


« S 


. 
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“Se IT IS FREQUENTLY, THE CASE THAT THE INVESTIGATOR WISHES TO 

. ASSEMBLE ITEM SCORES INTO CONTENT STRATA AND TO PROCESS THE * 

ITEMS WITHIN &ACH CONTENT- STRATUM AS ‘THOUGH THEY WERE ITEM 

POPULATIONS. PER SE’. CONSIDER» FOR EXAMPLE», THE CASE IN WHICH 

HALF THE ITEMS IN -EACH “SUBTEST MEASURE QUTCOME CASE A WITH THR 

OTHER HALF MEASURING OUTCOME CASE Be. HERE» THE INVESTIGATOR 

MAY WANT TO ANALYZE RESULTS SEPARATELY FOR A AND Be THIS 

° COULD BE. HANDLED THROUGH DATA STRUCTURES 1. AND 2. BUT TWO 
DUPLICATE SETS OF ITEM SCORES WOULD BE REQGIRED. TO AVOID 
USING DUREECATE SETS OF ATEM SCORESs- IT IS POSSIBLE WITHIN 

‘ THE PROGRAM.-T0' “STORE A DATA SET ON A SCRATCH FILE AND THEN. 

MAKE REPEATED” PASS ES THRUUGH ie laa OLE exert FORMAT CARULS® 

. (GF. THE © TYPE} ‘DESCRIBED IN 2C)e 


HASH ANAAARNAAAMANMAAFCANGAARSG 


e s 


THE ORGANIZATION OF INPUT CARDS FOR THIS CAPABILITY [S UNLY 
SLIGHTLY DIFFERENT FROM THE STANDARD CARD ORGANIZATIONse ASSUME, 
FOR EXAMPLE, THAT THERE ARE 5 PROBLEM SETS EACH OF WHICH IS TO - 
USE THE SAME DATA BASE. WITHIN THE DATA STREAM, THIS DATA BASE 
PRECEEDS THE PROBLEM SETS. THE DATA BASE MUST. BE PRECEEDED KY 

A CARD WITH "STR2" PUNCHED IN COLUMNS 1 THROUGH 4 AND TERMINATED 
BY A CARD WHICH HAS "END" PUNCHED IN COLUMNS ‘1 THRQUGH 3. 

EACH PROBLEM SET WHICH IS TO USE THIS DATA BASE MUST HAVE A, 
"1" PUNCHED IN COLUMN 14 UF CARD 2B. WHEN THIS UPTIUN IS USED, 
THE "DATA CARD SET" DESCRIBED IN ZC°JS OMITTED. é&VERYTHING 

ELSE STAYS THE SAME. . an 


Q 
io) 


A. PROBLEM SETS USING THIS STORAGE FEATURE: MAY. BE ae 
INTERSPERSED AMONG STANDARD PROBLEM SETS. 


- 


Be A CALL TO stone A DATA BASE MAY BE INSERTED AT ANY 
PLACE (PRECEEDING A PRUBLEM SET) WITHIN JHE DATA STREAM, 

A DATA BASE MAY BE STORED ANY NUMBER OF TIMES, 8UT 

« ONLY ONE BASE IS IN STORAGE AT ANY PARTICULAR TIMES 


OAPAMQAANNAANCAAANAO 


C. SCRATCH FILE 2 IS USED WITHIN THe PKOGRAM. THE 
STORAGE. CAPACITY OF THIS FILE 1S 13,650 CARD I MAUES. 
IF MORE STURAGE SPACE-IS NEEUED, AN ADDITIONAL FILE 


OOan 
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_ MUST BE-USED. — a . 
; m . 
C¥EKEK EXE hie ay HK eK eee reese ts *x iad JE EGOS IO eo x 


COMMON /BLOCKI/ SDATAt102+11)sNTEST 
DATA STR2/4HSTRZ/ : 
CIMENSION TITLE(18)° , 


READ spa RETR ECATIONS 3 FOR -PROBLEM SET 


1600. READ (5415£ND=5000) ° *(TITLE(1) 9121918) oi . 
IF ( TITLE(L) iNEq- STR2 ) GO TO 50 a 
CALL LOAD 
GO TO 1000:. 


50 " READ (557) K POP NTESTs NOLCHaNDIST AD ISK,UL,U25U3, us 
; IF (. NDISK .EQ. 1 ) REWIND 2 - -\.. 
., | WRITE (692) (TITLE() Ts 1418) ; 
ff IE { Ul, «LE. Of } WRITE (6,3) KPOP |NTEST »NOICHyNOIST 
ee S| ud eGTe Oo ) WRITE (694) KPOP sNOISTy U2 U2 93 U4. 
ki : 1 
CHECK ON PARAMETERS: es 4 a oe a 


‘NE RR=0- * 
IF ( NDICH wit. &. -OR. NDICH eGTs' 1 ) NERR=NERR#1- 3 
IF ( NOIST .LT. O° .ORe NDIST .GT% 3 ) NERR= NERR#1 
BE { KPOP .LE.°O ) NERR=NERR41 _ + 
IF ( Ul «NE. O. .AND. NIEST: .GT. 0 ) NERR=NERR#1 
aJF.( NERR oLE. 0 } GOTO 90 . 
' “ae ire (695) NERR ; a 
CALL EXIT, | % 
GQ CONTINUE = 


: BYPASS ESTIMATION PROCEDURE IF MOMENTS “ARE AVAILABLE ae 
(TE ( Ui .6T. 0. 9 GO-TO 100 
; ESTIMATE PARAMETERS FROM SUBTEST RESULTS 
: CALL. SUBTST( NDI CH KPOP »NDISK) 
C COMPUTE POOLED ESTIMATE eves Sane) OF4PARAMETER AND STANDARD 
C EFROR OF ESTIMATE. ASSOCIATED: Wi TH PARAMETER THROUGH JACKKNIFE PROCEDURE 
: CALL JENIFECULs U2 U3 USI 
‘ ESTIMATION OF RELIABILITY OF TOTAL TEST SCORE ON KPUP-ITEM TEST 
o . THETA= SDATA( 101, 8) : 
” RTT=KPUP*THETA/ (1. *k POP*THETA) 
| ) . WRITE (6,8) RIT - | 
; = GRADUATE HORNET AVE DISTRIBUTION ; 
) ‘ 7 13 


12 . = 2 . See 7 


° a wore 
at , 


100 IF’ 7 NDISt o£. 0 )° 66 TQ 1000" 
LE . -NDICH TEQe. 14) GO. -T0 1000 ce ; 
- S=Kpop . > ee < 
: cu TO ‘versus cee NDI ST s 
140° CALL NEGHGR {(KPOP,UryU2) .. - _* 4 4 
IF ¢ NDIST «LT. 3°) GO TO.1000 
120 CAEL PEARSH. (U1 9U29U37U4 900 9Sy1 0) 
GO. TO 1000 © 
-¢ eet . 4 
ik EXIT GRACEFULLY ers = 
re . 
$000 WRITE (626)- .° - i’ 
‘CALL EXIT ; ee 


Ne 


4 


c fica 
: FORMAT STATEMENTS te # is 


°F - FORMAT (18A4) : cag nhhy ; aa 

2 (FORMAT (1H1,18A4//) : 
-3 ‘FORMAT (31H SPECIFICATIONS FOR PROBLEM SET// - 

. 119H MAXIMUM TEST SCORET1I1/ oO 
. 219H NUMBER OF, SUBTESTSI11/- * : . 

‘313H ITEM SCORINGI1T/- a . - -_ 
4 ° ' G23H° NORMATIVE. ‘DISTRIBUTIONIT) :. a 
8 ty FORMAT (31H SPECIFICATIONS FOR PROBLEM SET// — - 

‘" VGH MAXIMUM TEST SCOREIL1/. : 
*223H NORMATIVE DISTREGUTEONIZ/ po. Ae ee 
34H UL"F15.5/3H U2F16.5/3H UZF16.5/3H U4FE6.5) — 2 

FORMAT (4H ***13,38H ERRORS ON CARD. SPECIFYING. PROBLEM SET//) 7 
- FORMAT (20H ALL INPUT PROCESSED) sd 

‘FORMAT, (2145 31291X,4F15.0) “ 

END (7/454 ESTIMATED RELIABILITY | OF TOTAL TEST SCORE ISF 8% 5/1) 

EN a. : 


SAYA 


. 
{y 


‘¢ dl 


"tes 


@ - SUBROUTINE LOAD te . F 
ae _— G 4% me 
"DUPLICATION or DATA ‘BASE “ON SCRATCH FILE. NUMBER, a 4 


—_ DIMENSION. C200 ts  # oe 6 
; - DATA END/ 3HEND/:.« .. aan ; ge 3 
"REWIND. 2 - ee ae ‘se 
—— WRITE. (642) 7. eo. se." * 
. 100 READ (S411: (XU), Ted ,208; ewe 2, . 
n WRITE <6, 3) {KUL} p1=1920 at : ‘ .:, 
“°F @ XL) .EQ. END-3'GO TO 200 ": . : 
"WRITE (291) (XOD) st= 1920) é- * 7 oe 
GO T0100 
200 END FILE 2 ; . 
RE TORN ; - 4 . * 
1 ° FORMAT (204) 
. 2 FORMAT (//37-HICARD IMAGES LOADED ON SCRATCH FILEs 2//) 
f ge oe FORMAT | 1X2 044) og : — ba 
pete END ° ’ a i. a 


Pa 


z?- 


7 
. 
ae, 


. @ SUBROUTINE NEGHGR (POP, XBAR »XVAR) pes — , 4 


C ay 
C _ GRADUATION OF NORMATIVE DISTRIBUTION BY NEGATIVE HYPERGEOMETRIC as 
C. OF STRIBUTION a . mee ‘ , 5 
c -," os . Ree . : ad 
oC. cb MPUTATION' OF CONST ANTS = : an Be a 
C 3 * 
a6 A21=4KPOP/ (K POP ~ te 140 14 =xBAie( Kop xo aRD /tKPOPHAVARD e 
WRITE 4691) XBARsXVAR, A2l 
‘2 . SF. ts “A21 GT. 0. eo ANDs one rath 1e ) GO TO 120 ‘ 
| WRITE (6920 ae 
* RETURN’ |; - "* - 
: 120 A=(-1.41. /A21)*XBAR a ‘ & 4 
2 Tet, Ban Le +KPOP/A21 a ssl 
7 ".$LOG1=0. = oe a . 
/ « S$LOGZ=0. , : ; ; : - ‘s 
‘ oy CAB 2 ae 
3 00 140. 1=1 sk POP : 
- SLOG1= =SLOGI+ALOG10(8=141.) ee op 
140 SL0G2=SLOG2+ALOG10(C~I+41.) - i J a 
C= 10.##( SLOG 1-SLOG2) & sa .: . 
c — : . 
C GENERATION OF NEGATIVE. HY PERGEOMETRIC D IS TRIBUTION 


. 


* u 


WRITE (673) eo i x P 
S=KPOP 


-@ > N3SKPOP+1 . 


©, 


1 


GO TO 165 
HX=C 
Piaaaniibiles 
“CONTINUE 


"RETURN * 


END 


wt 


-00 Oass=2)N3 


DG /L50 I= K--... 
_SLOGI=S£ 06410 LOS-141.0) 
SLOG2=SLOG 2+ ALOG10(A¥1-1.Y 
SLOG3=SLOG3+ ALOG1O{B-141.) ; ~ 
SLOG4=SLOG4t+ ALOGIO(FLOATCI)) 
“HX=C®10. #*(SLOGI+SLOG2- 7SLOG3-SL064) 


WRITE (694). /KoHK 4K 


. 
“~~ . 
? 


K “EQ. 9. My Go" to, 160° 


‘*. 


FORMAT (//82H APPROXIMATION OF NORMATIVE DISTRIBUTION WITH NEGATIV 

1E, HYPERGEOMETRIC, DISTRIBUTION//21H ESTEMATED parr ener en MEAN,» | 

*» 2F25¢5/9H VARIANCEF 21 .5/4H A21F26.5) - - 
FORMAT... (//30H, DISTRIBUTION ABORTED Adil LE O//) 

FORMAT (175% 5HINDEX #5X91 3HRELAT IVE FREQs 10X12HCUM REL FREQ//) 


FORMAT {1XyI9 a1 85, F22.5). 


a 


"SUBROUTINE TABLE (Ts Ny AREA) . ; ' 


CALCULATION OF STUDENT®S T CORRESPONDING ‘TO PRESCRIBED ALPHA aetee 
(TWO-TAILED PROBABILITY} 


ARGUMENTS 
- T. ~ 
N 

AREA 


STUDENTtS T 
DEGREES UF: FREEDOM ASSOCIATED: WITH 1. « 
RELATIVE AREA aca ha WITH CONFIDENCE INTERVAL 


hou the 


POLYNOMIAL APPROXIMATIONS USED FOR ‘BOTH T ‘AND. Z 


‘REFERENCE: ABRAMOWITZ, Me. & STEGUM, I. A. HANDBOOK. OF MATHEMATICAL 
‘..> FUNCTIONS. NATIONAL BUREAU OF STAND ARDS APPLIED 
.'” " MATHEMATICS SERIES NOe 55, 1964. ’ 
* (SEE 2622.23 AND 26.7.5) ° 


- TF “CON 2G6T. 0 2°60 TO 100° 0 i 
WRITE, (691) N 2 ook aes bg 1 Re 
GO To 150 * ae ots 

100 “P=(1e-AREAD/22>.... a4 3 a te 

IF { P GT. 0. - AND. Poles 5) GO TO. 200 i oe 

WRITE (692) P. . ry a 


i 150 T= O. " f : ae ; Ts, 

r RE TURN h eo ee, : 
co X=SQRT(ALOGE 2e/ (PHP)? ee ty 

z= X- (26 51551.74. 802 853% X+ 6 O0L0328#X#X) / ‘ 
l gtboeteg7eTBOhar ol Speco" At 00] 308% xe#3) 

L2=1*83 et : “~ ‘wi a. 
L5=LT#*S 7 7 
(7=(**T 
LO=LEEO, ; : 
GL=(Z342)/40 * > 4 a [= 
G2=(5.*Z5416. #13 435*1) 96s ! oe 

G3=( 3.427419 o#Z 54170 4Z 3-156%2) /384. 

G4 =(T9e#L9*T 160%27 +148 20 ¥Z5¢ 1920042 3- -945. #1) 192160. 


T=2+G1/N+G2/ N*#Z+G3/N#*34G4/NOH4 ce : 
RETURN a re — 
FORMAT (//27H ERROR tN T~TABLE ARG: OF =112//) - fy : 
“FORMAT (//26H ‘ERROR IN: I= a ARG: P =E15.5//)9 

END ‘ 


17 


° ca > a 16 : ~ ai 
"SUBROUTINE ROFMT (FMT2N) 0 a * #: 
SUBROUTINE FOR INPUTTING VARIABLE. FORMAT a 
INPUT STRUCTURE »— be te — a | 


FORMAT (ENCLOSED IN PARENTHESES) - COLUMNS 1- 72 
CONTINUE ON-CARD 2 IF NECESSARY 
“CONTINUE ON CARD 3 IF NECESSARY 7 i 
CONTINUE TO CARD 9 (PAX) IF NECESSARY 

* “END cas FORMAT" RUNCHED IN a 13 


"AANMAMAAANAADQO: 


DIMENSION. FMT (162) “ , F 
DATA’ ENDyBLK/3HEND y1H / eS 
DO 50 1=1,162 
50 °FMT(1)=BLK -_ 
“N=l i. . . ‘ 
DO 100 I=1,10° ee 
* MENtl7 : 
-, -RRAD (5,1) CENT) yd =NyM) 
.  , IF (:FMT(N} «EQ. END-) RETURN 
100 .N=N+16. | a re : tie 
WRITE (692). - 2 got ee ee 
STOP 2 eee 
FGRMAT (18A4) ° 
FORMAT (37H 6 EXCESSIVE NUMBER OF FORMAT CARDS} ae : 
END ~ ae ol is 6 


"OND bee 


FUNCTION FRV (RsVV) oo _ 4 


EVALUATION oF COS (X)*#R EXP (- V#KIOX FROM “PI fe TO PIs2 USING 
SIMPSON? S RULE 


CV OY > 


V=ABS(VV) ‘s 4 

PI=3.141593 , 

O=PI/1G0.  * te 

FRV=06 ‘ .s . . . , ZS 

DO 100 12149942 ° “4 YS 

S-1.570796414D te 

HRV=FRVt2. #COS(T )AHRSEXP (VAT) #COS(T 45 1 #RHE XP (> VT=VD) 

ER V=2.*D¥PRV/3q 

RETURN . ; 7, - & | a 
END ; , ee .. 4 ; * 


/ : F \ . _ saan ; ; . ; 


POLYNOMTAL/ APPROXIMATION 


“Maan 


; TF. -{ X (87, } 62694 _ 
4 WRITE ((by2}? - | 


GAMMA#1.E75 .- ; 
RE TURN e 
1 6 K=XX/o oa 
ie 2E-6 
; rae X26) 50950915" . 
10 TF i{ X-2. 210,110,915 A 
15° X=Xe1e . Z 
© GHEGKEX 
601010 | a . 
50° IF (-X=le } "60+4205 110/ i 
G 
C . Test FOR, -X NEAR NEGATIVE INTEGER OR ZERO ; 
ce : 
60 jiF ( XRERR.) 62, 162480 ze : : 
GZ: | Y=FLOATCINT(X))-X- 
| IF., / ABSCYI-ERR*) 1304130564 | 
et] FRC Le-¥-ERR i 130) 1304 70 a 4 
Cc . ! 7 
x Not NEAR A, NEGATIVE INTEGER oR. ZERO “ 


€ 4 wi 
ae iF ( X-1le b 801804110 : 
6&0 GX=GX/X- 
X=Xtle 
- GQ“TO 70 
740 Yex~ 1. : 


* 


1{-.56847294+V4#(. 2548205+ 8 {-20514 930004000) 
GAMMA=GX4GY, 
120 RETURN . 
130 WRITE (691)5 
RETURN ° 
L FORMAT (/53H GAMMA ARG WITHIN: -00 001. OF. BEING. A NEGATIVE INTEGERS 
- 0) 
Zz FORMAT (/44H GAMMA ARG GT 57, OVERFLOW, ANS. SET TO 1eE75/) 
-END a * a7 is 


s° = Pa . : ; 18 . : : rs : a af 
= ‘ * ‘ . : | . 8 ; ° ‘ 
SUBROUTINE JKNIFE (U1,U2,U3,U4) * 
‘ ' a 4 
CALCULATION OF POOLED ESTIMATE GF PARAMETEP ANU ASSOCIATED STANDARD 
ERROR OF ESTIMATE (BY THE JACKKNIFE PROCEDURE 


MAP FOR SDATA MATRIX 23 4 


SDATA (491) TC SDATA(J 8) , PARAMETER ESTIMATES FROM SUBTEST J 
SDATA(J,9) | NO. OF} EXAMINEES FOR SUBTEST J 
SDATA(J, 10) NO. OF|ITEMS FOR SUBTEST J 
SDATA(J511) ‘ + NOs OF OBSERVATIONS FOR SUBTEST J 
Sa ea POOLED ESTIMATE FOR PARAMETER. I. 
SDATAX102, 1) 4 ee ESTIMATED: SE FOR PARAMETER I 


i 
= 
‘ ‘ 


GOMMON /BLOCKI/ SDATAC 16211) NTEST 

DIMENSION A(8)_ B(8)y PRCNT (6) sCONFE 6) 4 

DATA CONF 4H70.014H8060, 4H90.0% 4H95.00 5 4H9T« 544H99. 0/ 
DATA PRCNT/e 70, 80; ge 90 ge 95 90:575 90997 ¢ " 


MAAADANAANAAMO 


IF { NTEST eLEZ 1 } RETURN 


DD 200 I=1s6 a : 


€ 

c ae. 

C CALCULATION OF WEIGHTED PARAMETER ESTIMATE USING-~ALL SUBTEST RESULTS 
* C ‘ ‘6 2°” : 


PALL=0. 

.TOBS=0. . 

DO 100 J=1LsNTEST a: 

SDATAC Ss 11 F=SDATAC Js 9) *SDATA{J,10) 

TOBS= TGOBS+#SDATA{ Js1l) : 
100 PALL= PALL#SDATA( Jy I)ESDATA( Jel) : 


~ 


C 
C CALCULATION, OF PSEUDOVALUES PER SUBTEST Por A GIVEN PARAMETER . 
C ae 
‘SP =0. - ' 2 7 
- SPP=05. — : a 
DO 150° J=1yNTEST - ‘ 
P= =NTEST#PALL/TOB S- (NTEST#1 «) #( PALLr SDATAC Sy] PRSDATAC Sy 1100/ 
1 TORS Pav Atset 2), 
- - SP=SP+P .° ~ j 4 - j aS 
150 SP P= SPP+P*P He ‘ _— ‘ 
C - . 


C CALCULATION. OF POOLED ESTIMATE OF: PARAMETER AND ESTIMATE OF STANDARD 
C ERROR OF ESTIMATE ASSOCIATED WITH RARAMETEK 


C 
SDATA(101 I) = =SP/NTEST 
SDATA(10291)=0. 
T= (SPP-SP*SP /NTEST )/ (NTEST*(NTEST~1.)) | fs iy 
IF ( T .GTs 0. ) SDATA(LOZs1)=SQRTIT) 


C , . ioe 


Lt — 
200 CONTINUE 


¢ 
Ct, 


* ah * * oe ¢ - oto i “4 ete ey 8 4d es 
ca : ; *, 3 
.- 19 - : . 
fo e = 7 at 
* : ‘ 
ae a’ . 
’ 


a“ 
_— 


WRITE (Os1¥ ((SDATAC Jy E) g1=198)yJ=101,102) 


: ; oa 
-CGOMPUTE CONFIDENCE INTERVALS FOR PARAMETERS 
(DF=NUMBER OF -SUBTESTS HINOS ONE) 


Nana - 


DO 300 J=3,6 
CALL TTABLE (T)NTEST- LyPRONTIJ)) 
WRITE 4692) CONF (J) 
DG 250 1=1,8 
TEMP=T#SDATAULO2Zy 1) 
ACI) =SDATA(10191)-TEMP oe 
250 -BCI)=SDATA(101,1)+TEMP t- 
WRITE (693) CACK) 9K=118) 9 (B(K) sK=198) 
300 CONTINUE ~ 


UL =SDATA(101 91) 
U2=SDATA(101 92) - 
U3=SDATAC 101 » 3) : /& 
US=SDATA(1O1 +4) : | —_ 
i . ve , : . ; . ” 
i - RETURN ‘ - 
Cc " th. : 
1 FORMAT (//5H MEANGE15.5/5H SE B8E15.5//) ~ 
gy |... FORMAT {1%sA4y1Xy27HPERCENT. CONFIDENCE INTERVAL) 
@ 3. FORMAT (5H FROMBELS.5/5H TO 8615.5//}” 
3 END | 


SUBROUTINE PEARSN(UL sU21U3,U4sL19L2y STEP) 


x a 


“ANNAN MEOAMANANAAAMABAAAA.A 
c 
~~ 


GRADUATION OF F REQUENCY DISTRIBUTIGNS WITH PEARSON CURVES 
PRIMARY REFERENCE 


ELDERTON». We P. AND. JOHNSON, NeLe SYSTEMS-OF FREQUENCY CURVES. 
“LONDON/ CAMBRIDGE UNIVERSITY PRESS/1969. (PAGES 35-109) 


’ 


e 


ARGUMENTS OF SUBROUTINE 


. MEAN 
U2 2ND CENTRAL MOMENT (VARIANCE) 
U3 3RD CENTRAL MOMENT 
U4 4TH CENTRAL MOMENT 
tl LOWER LIMI NDEX 
L2 UPPER LIMIT OF INDEX 
STEP DESIRED STEP SIZE IN GGING FROM .L1 TO L2 


‘ 


REQUIRED. SUBROUTINES 


~~ 


GAMMA -EXP .ABS SQRT* COS ATAN 


REAL MyM1l M2 et1,bk2,L AMBDA. : 
DI MENS ION CELA RRL O TRL 


Cc 
Cc COMPUTATION OF INDEX OF SKEWNESS teL) AND DEGREE OF. KURTOSIS (B2) 
C AND CRITERLON VALUE ‘USEO on THE SELECTION OF APPROPRIATE PEARSON 
C CURVE _ 
C : Me 

X0=0. 

XN=0. 

NHIGH=0 

CUMP=0. 

P1=3.1415927 < 

PPI=2.¥*PI : 

Bl =U3*U3/(U2**3) 

as B2=U4/ (U2*U2 ) : ; 

Z1=B1 . , ‘ ‘ 

Z2=B2 

CRITI=(B1*(B243 . )¥*2 )/(4.4(4.*B2- a *B1)*(2.*B2-3.*B1- 60)? 

WRITE (6914) Ul.U29U3,U4¢ 

#RITE (6,91) B1l,B2,CRIT1L 

If’ ( U2 GT. O- ) GO TO 60 

WRITE (6,15) : 

RETURN ens ee 

60 IF ( Li -LT.~ L2 «AND. STEP GT. O- ) GO TO 65 
WRITE (6,16), : 
»RETURN 
65 CUNT IENUE 

A=4.*BZ-3-%*B1 . 
@ B=3.*B1~20*B246. 

CRIT 2=(A*(10.*B2-12. *B l- 18.) ##2- B1l*(B2+3. Ieszate. *B2-95*B1- ~“l2.))/ 
1 (B¥(BL*(6243. )**24+4, ¥A%B)) 


eC 


ered 


a 


e@ i 


aan 


roe kay 


AON 


21 


CRIT 3=5.¥*B2-6.4B1-9. 
CRIT4=2.*B2—3.*Bl 66 
WRITE €6)23) CRIT2,CRIT3gCRIT4 


TEST FOR VALUES OF CRITERIA APPROXIMATING THEORETICAL VALUES 
-, IF ( CRITLLE. 202 «ANDO. CRITI.GE. ~.02 ) CRITL=0. — 
IF € CRITL-LE. 1.02 .AND. CRIT1-GE. -98 ) CRITI=1. 
IF ( CRITZ sGE-e -.02 .AND. CRIT3 eLE. 202 } CRIT3=0. 
IF ( CRIT4 .GT. ~202 eAND- CRIT4 oLT. 202 ) CRIT4=0. 
IF ( Bl -GEe 3498 AND. 81 eLEe 4-02 ) 21=4. 
IF ( B2 Le 3-02 eANDe B2 -GEe 2-98 ) Z22=3. 
TF ( B2 -GEe 8-98 eAND. BZ eLE. 9-02 ) 22=9. 
TEST FOR CURVE BEING MAIN TYPE GR TRANSITION TYPE 
AF ( CRIT4 «EQ. 0. ) G6 TO 600 - 
IF ( Zl oEQ. 4. ‘a AND « zZ2 eEQ. 9. ) GO TO 800 
IF ( CRIT2 ) 8097980 , : 
‘79 IF ( CRITL eLT.e O« e AND. CRIT3 wLT. 0. } GO TO 700 
IF ( CRIT1 -GT. 1. «AND. CRIT4 .GT. 0. ) GO TO 850 
IF ( CRITI.LT.O.eANDCRIT3.GT. DeeAnh eC LT Ys LT.0- ) GO TO 750 
80 IF ( CRIFL) 1009400,90 ~ : ; 
30 IF ( CRITI-1.) 200,500,300 ’ 
CURVE I *FIRST MAIN TYPIE* 


i00 WRITE (6,4) . 
R=6.*(B2-B1—1.)/(60+3-*B1-22*82) 


C= eee 1*#24 166 *( REL ) é 
AL2=SQR Ci/2. - 

C=R*(R+2 TAS ORT (51/0) \ 7 
M1=(R=-2e-C)/2. . 

M2=(R=2e4CV/2. . 

IF. U3 IS POSITIVE, M2 1S POSITIVE ROOT 

IF ( U3 eLE. Of ) GO TO 120 

IF ( M2.eGEe Ml } GOTO 120. b % - 
TEMP=M2 a on 
M2=M1 : oo 
M1=TEMP . 

120 AL=(M141.)*A12/(M14M242.). 

A2=A12-Al - ° 
VE=(1e/A12 EC ( (M1410) HEM LEC M2416) M2) /(MLEM242.0 ) (M1 MZ) 
1(GAMMA(M14+M2 42. )/(GAMMAC M141 .)*GAMMA(M2#16))) — - 
WRITE (692) ; 

X=Lh : ‘ 

KK =] . : _ 

150 IF { X GT. L2 )} GO TO 1000 
T=X-U1. ; : eo 4 
R1=1e+T/AL : 


R2=Le~T/A2 


~ 


4 . IF ( Ri eLE. Oe eAND. R2 oGTe Oe ) »P=VYEeRQ*K M2 
7 © IF ( R2 slEe Oe eANOSRIL .GT. 01 2) P=YE*R1¥*M1 
‘ . TF ¢ Rl. eLE. 0. eANBs. R2' eLEe Of1 +) P=YE 
IF ( R1 SGT. Os eANDe R2°2GTe Oe ) P=YVE*R1] €*€ML¥RZKKMD 
“IF ¢ R11 WLE. Oe eORe R2 eLE. O01 } NHIGH=1 
GO TO 5000 , 


Cc 
C CURVE IV *SECOND MAIN TYPE* 
s¢ ? _ 
200 WRITE (6,5) 

R=60*(B2-Bl-1.)/ (2.* 82-3 .#B1-6. 

M=(Rt2.)/26 
: , A= SORT((U2/160) #( 160 #( R= ieycauetesgoieesis: 
ve SQRT (URER* (R-2.6 )4**#24B1)/(1604(R-1. )-BL¥*(R-2.)*#2)) 


C 
C V AND-A HAVE OPPOSITE SIGN OF U3 
Cc ; ‘ ~ 
SGN= 1. ~ 
IF ( U3 eGTe O« ) SGN=-1. 
A=SGN*ABS(A) , . : 
V= SGN*ABS(V) a : 
C . 
C NOTE see ARCTAN (-Z)=-ARCTAN(Z) 
C . 
THETA=ATAN(ABS(V/R)) 
TF ( V/R eh Te Oo ) THETA=~THETA 
COST=COSCABS (THETA)) 
: 
oo . . 
IF ( V-2e ) 225,225,250 
ae ; 
225 NOs bod (ROSCA 04ERVORS VD) 
' GO-TO 270 
C a. 
C APPROXIMATION FOR YO. 
t ; . ’ : 
250 FEMP={COST , — ##2)/(34#R)—1e/(12.*RI-THETA*V : 
OSC /ADASQRTUR /TPEIHEXP (TEMP) /( COST _ 8#(R#1LS)) 
C <a 
C ee. bi % ; : ; “ 
270 CONTINUE “ : 7 
eee Gee) . 
=Li 
YEAR 
om _KK= =2 , 
/ 80 LR CX .6T. 2. } GO TO 1000 © 
. " T2X- U1 a , Z 
- 7. Sty A-Y 
So PzVOe(1. +0084 (- M)REXPC- ~V#ATAN(Q)) 
GC TO, (pooo- : 


C 
e: CURVE vi “#THIRD, MAIN TYPE* 


7 iii . 
: pe - 


PERG: © ge * # i a4 


"300 WRITE (6:6) - ‘ 
Q] — .--Rr6.#(82-B1- 1.9/1 6.+3.481-26482) e 
; ~ * CHBL¥(RH+2 2 )##2416.4(RE1.) : 3 
. A= SQRT(U2*C) /2. . is 


A HAS SAME SIGN AS U3 


OOM 


IF ( U3 eLT. Oe } AS-A 
TE MP=R¥ (R42. )/2.*SQRT(BL/C) . 
Ql=-(R=-2.)/2.+TEMP 


Q2=(R-2.9/2.#TEMP : 


aan 


_R IS ALWAYS NEGATIVE HERE. Q1 IS GREATER THAN Q2 ae 
. TL=Q2+t1. 4s ° in. 
LL T2=Q 1-Q 2-2 e 2 
fe T3=Q1l-1. 
° AL =A*T3/T2 - 
A2=A4T1/T2, ee ’. 
YE=(T1¥¥Q24T 24% (Q1-Q2) #6 MA(QL))/(ABS(A) 
1 +#T34**Q1*GAMMA(T241 « )*G AMMA (T1)) 
WRITE (692) _ : 
X=L1 7 ae oe Ke 
KK=3. , ' es 
350 IF ( X .GT.e L2 ) GO TO.1000 _ = =? » %e 
_ ‘T=X-Ul ae 
© R1=ABS(.1.+T/Al) \ 
RZ=ABS(1et#T/A2) . 
P= VERL**I-Q1)#R2¢*02- 
vO TO 5000 


TRANSITION CURVES 
400 IF -( Z2 -36 } 4501410,480 ' : 
NURMAL CURVE OF ERROR *TRANSITION CURVE* 


410 ‘WRITE (697) , - 
YO=1 o/SQRT(TPI*U2) | - " a 

WRITE (652). ) es 

X=L1 

KK =4. ° : . are 

420 ate ( X .GT. L2 ) GO TO 1000 

: = YOREXP (—( ( X-U1)#*29/ (2 6*U2))- 

60 TO 5000 


ae ‘ 
* € TYPE II CURVE * TRANSITION CURVE* 
+ C = , 


~ 


' 450 WRITE (6, 8) 


M=(5.*B2-9. 1/(60-20*B2) s é 
e « R= (2 640282 11 (3 «82 A 
@ . =(1./SQRT A2¥P1))*GAMMA(M41. 51/GAMMA (M41 ). 
an ITE (632) Yi 


¥ " _ Kz LI : , = . * 


KK=5 : 
O IF-{ X .GT. £2 2 GE TO 1000 
_ T=X-U1 , 
RL=1 e-T#T/A2 .- 
IF ( RL «LE. 0.:) P=YO 
IF ( Ri GTe. Oo ) P=YO*R1**M 
' 60 TO 5000; 
Cc : 
C TYPE VII CURVE *2RAnN 
C 


ye’ WRETE (6,2) 
Be X=11 ‘ . ¢ 
450° IF { \. GT. L2 ) GO TO 1000, 

T=X-U1l 

P=YO*(1 4 T#T /AZ) ** (-M) 

GOTO 5000 


GAMMA(M)/GAMMA(M- 25) 


TYPE Vv CURVE *TRANSITION con * 


"500 - WRITE (6,107 
: P=4ot(B.t4.#SQORT(464BIV 1/81 _ 
. GS(P=2.)*SORT(U2*( P-3Z)) 
YO=G¥#*(P—1.) /GAMMA ( P- 
ORIGIN=U1=G/ (P-2.} ; 
“TF AC U3 olTe OY ) ORIGIN=U1FG/(P-2.) 
WRITE (692) 
-X=L1 
KK=7 
IF (.X% eG6T. L2 ) 0 1000 
X= X-ORIGIN 
T=ABS(X) 
P= YOX TH (—P) KEXP (-G/X) 
GO TO. 5009, 


4 


ae 


TYPE ILI CURVE %TRANSIT ION CURVZ* 


600° IF ( Bl sEQ. 0. ) GO TO 80 
/ WRITE (6511) 
LAMBDA=2.*U2/U3 
F=4./Bl-1» 
A=(F+1.) /LAMBDA 
YE= LAMBDA#(( (F41.)*#F) /(EXPUFEL« DEGAMMACE4+16))) 
WRITE (692) ; 
X=L1 | 
KK=8- 
If ( X .6Te L2 ) GO TG 1000 
T= XU] : 
KL=1e4+T/A 


IF € RL eLEs oF ) p= > YEXEXP(~ LAMBDAST) 
AF ( RL oGTS°0. ) PSYESR14*F SEXP - LAMBDAYT), 
| TF ( R1 “LE 06 j NHIGH= i 
“G0 To 5000. f/f td 2 5 
TYPE VII CURVE *TRANSETI ON CURVES eS 4 ; : 
“700 WRITE (6519) 
_ SOLUTION. OF CUBIC EQUATION FOR VALUE. OF M 


WITH TYPE vill, M IS RESTRICTED TO 0-3 RANGE 
“CF(L) =160*B1 a ey oe © ee 
- CF(2)=524.¥*61 =: ae ‘ - 
" . CF(3)=96*B1-125 : is 
CF(4)=4.-B1 | : 
CALL POLRT (CFr39RRyR1y TERY -y 
DO 705 I=1,3 . ‘ 
M=RR(T) 
~ IF (MeGEe0eeAND. MeLE +12 ANDSRI(T) +6004 ) GO TO #10. 
705 CONTINUE to 
WRITE (6525) te mo ae 
GG. TO 750 Ps a ae 
“710, ie (2-H) #SQRT(U2" (3. “MI /(1.-H) 
A=T#(Le-MI/E20-M) a 


eae oe. 


AO 


A HAS OPPOSITE SIGN OF U3 - 


IF { U3" elTe 0 ) A=-A 

YE=C( Le MD*( 20- nah 

WRITE (692) 3 , 

2. UO XSLD - eo re a 
: KK=9 oe ede oe ‘ 4 
‘720 .1F { X eGTe L2-) GO TO 1000 - . 

 TeXeUL 

P=YE#(14—T/A)##(-M) 

GO TO 5000 _ , * P 


C TYPE IX CURVE *TRANSE TION: CURVE* 


750. WRITE (6920) t a 


) 
& 


. SOLUTION oF CUBIC EQUATION FOR: VALUE ll M 


WITH TYPE. 1% M MUST BE NON- NEGATIVE 


anal, 


. CF(1)=160%B1 
CF(2)=24-*B1 A 
CF (3)=9.¥B1-12, 
CF (4 )=B61-4. 
. CALLE POLRT (CFp3yRRyRIZTERD 
~ 2 00 755 T=l93 
| MERC 


vagy “> 
vs 


e 1 
Pe 


a 2S a ae ar . ~ 
oS grr os! : 
@ | ae 
a 
er 
x 


AOBAaAA lacga ~ 


vw ~ 


Cc 


800 


820. 


850 


WRITE (6926) 

60 TO “100 P23 

T=(Mt26 eS ORTUZ# (MES. 17 iMeD vy 
- A=(MtI. Y*T/OMt2.) 


YE=(M41. D*FCMFL IAT HUME 2. em) 


WRITE (6,52) 

X=t1 - 
KK=10 * a 

{F ( X.-GT... L2' ) GO ‘TO 1000 - 
T=X-Ul ; 


P=YEX(1.4T/ADHRM 


GU. TO 5000: .. 


ee TYPE. X CURVE *TRANSITION CURVES 


WRITE (921) 


> A=SQRT (U2') 


YE=1./(2.718282*A) - 

WRITE (632) 

X=L1 : 

KK=11 

IF ( X .GT. Le ),GO- To 1000.. 

T= (X-u1) 

IF ('U3 LT. OQ. } P=YE*EXP(T/AY 
“LF ( U3 .GE. QO. ) P2YE#EXP (~ Tsay 


GO 0 5000 - coos, 


TYPE. xi CURVE TRANSITION curves 


WRITE’ (6,22) 


rd 


wiTH TYRE/ KI'y M MUST BE. GREATER THAN 1 


° 


CECL)= =16.*B1 

CF (2)4-24.*B1 

‘CF (3)=9.¥*B1-12. . 
CFL4)=4.~-B1 oe * 
CALL POLRT (CF 3yRRe RTs TER) 
DO 855 \I=1 43 

M=RR (EF)! . 


ye: ¢ 6B. Os sANDe- RICE) EQ. oe) 60 
3 CONTINUE ©. 


RF { M .GTe le “AND. RIT) EOS Oe. ) 


- .CONTENU ba y’ 


WRITE (6424) 


60 TO 0 
* B=(M-2.) ASQRT(UZ*(M-3- / (M16) } 


YE=(€M~2.)**M)/(B¥(M—-1.)*#(M-1.)) 


HAS SAME SIGN AS U3 


B=ABS(B) 


e 


‘SOLUTION OF. CUBIC EQUATION FOR VALUE OF M 


“GO TO 869 


eee G8 Ey Oa ee is at a eS 
ae z get ge Wie ots eee 
ae + - Oe "tf Mgt 
me + a ast . 
= ac . - = sf 
ere oe ee = 
or = . tee. 
. € ee jg . A 
s . bs * 
. wo, > . 
. +; * 
> » te ‘ 
T > 7 i: 
. . - 
‘ Bn de. s 
fa) 
4% 
4 ~: 
4 
7 t 
° 


TF { UB «LT. 0. ) Be- 
A=BECM-1.)/(M=20) 
WRITE {652} 


X=L1 . : Sw : me ge si 
; KK=12 : : s @ * oe a ; aa 

8-70 IF ( X .GT.e L2 } GO a 1000 _ _— Se , > 
’ T=X-U1. : “s . ‘ e, Rape 4 

P=YE*(1. +T /A)**(- -M) - - a : 

GU TO 3000 : : ae id a : 
C "ad, 7 = oe a 
C AREA UNDER’ CURVE SEE ROMINATER BY TRAPEZOID FORMULA . 4 4. 


C os. 

1000 Are # (CUMP~(XO#XND /2 a, a a ee 

CARTE (6913) AREA | |, ar a ; 
RETURN. es e & a2 


oe a 7 a, 
Cc 2 ‘ ‘ fr -y ‘. 
* 5000 WRITE (693) 'XoP a 7 ‘ 
IF- ( P .GEe le ) side ~ F 
.  CUMP=CUMP+P 2 bg s & . < 
- IF ( NHIGH .EQ. 0, 560 To 5010... - oF yO #2. 
WRITE (6,18) a , a a 
‘NH IGH=0 a a oe re : oa 4 
“ —« 5OLO IF €( X 2EQe LL ) AOR et MS es > ee. o 
~ - (TF € X EQ. b2) XN=P oe —_ - ie 
X=X+STEP ike ’ 
GO TO (1501201550420 460,480,520) 6204720,726)820/870), KK - 
eo © 2. 4 & - 
C : - . 
1 -FORMAT €° / 16H" INDEX OF SKEWNESS F2 165/194 Decreé oF KURTOSIS- 
'  1F20e5/16H CRITERION VALUEF23.5)) -  -; ; =. 4 nm 


2 FORMAT (5X,5 BINDEX » 5X; A3HRELAT IVE FREQ//) 
3 FORMAT '(1X9F9.35F18.5) ". cet? 7 
+4 FORMAT (/22H TYPE 1 CURVE SELEGTED/) . > Pa cee a 

5 FORMAT (/23H TYPE IV CURVEsSELECTED/) 
,FORMAT (/23H' TYPE VI CURVE SELECTED/) 
7 FORMAT (/28H NORMAL ERROR CURVE SELECTED/)- -.,° 4 
8 FURMAT ('/23H TYPE LI CURVE SELECTED/) 
9 FORMAT (/24H TYPE VII CURVE SELECTED/) = \ ; 
16 FORMAT (£22H TYPE V CURVE SELECTED/) 
11° FORMAT (/24H TYPE JIT CURVE SELECTED/) 2 
. 12 FORMAT (//15H APPROX AREA ISF15.7//) © to 
"14 FORMAT, (////49H APPROXIMATION GF DISTRIBUTION WITH PEARSON CURVE// 
.  12iH ESTIMATED PARAMETERS///5H MEANF34.5/9H VARIANCEF3005/19H 3RD C 
' 2ENTRAL MOMENT F20.5/19H 4TH CENTRAL MOMENTF20.5) , 
15 FGRMAT (//23H NO CURVE WHEN U(2) = O//) ~~ ; 
16 FORMAT (//38H INVALID LIMITS.AND/OR DELTA ARGUMENTS//) 
18 FORMAT (1LH#,30X,32HORDEINATE UNDEFINED AT THIS POINT) 
19 FORMAT (/25H TYPE VIII, CURVE SELECTED/) 
20 FORMAT (/23H TYPE IX CURVE SELECTED/) 
21 FORMAT (/22H TYPE X CURVE SELECTED/) s 
22 FORMAT (/23H TYPE XI CURVE SELECTEO/) — i 
23. FORMAT (/7H LAMBDA, E36.5/10H 562~6B.-99F2925/10H’ 2BZ-3B1—69F2965/) 
é4  FGRMAT (/22H M,eLte le FOR TYPE XI /) } 


28 — ~ 


7 wo 25° FORMAT (/29H M QUT OF RANGE FOR TYPE VIII fy “_. = 
/W 26 = FORMAT (/23H M NEGATIVE FOR TYPE IX /}. - 2 * ee 
END 


; “s ; SUBROUT INE_POLRT (XCOFsM4RODTR;ROOTI y LER) 
; COMPUTATION OF THE REAL AND COMPLEX ROOTS oF A REAL POLYNOMIAL , 


aan 


METHOD? NEWTON-RAPHSON ‘ITERATIVE. scchnigle., “THE FINAL ITERATIONS ON 
EACH ROOT ARE PERFORMED USING THE ORIGINAL POLYNOMIAL RATHER 
THAN THE REDUCED. POL YNOMIAL TO AVOID PECUNUER TED ERRORS IN THE 
POLYNOMIAL 


a7 


LEMITAT IONS: LIMITED TO 36TH ORDER POLYNOMIAL OR LESS. 

Ps FLOATING POINT OVERFLOW MAY OCCUR FOR HIGH ORDER 
POLYNOMIALS BUT WILL NOT AFFECT ERE ACCURACY OF THE 
RESULTS.« 


PARAMETERS : o$ 


A tore, an 


XC OF VECTOR OF +1 COEFFICIENTS - OF THE POLYNOMIAL ORDERED FROM 
a, SMALLEST TO LARGEST “POWER 

COF | WORKING VECTOR OF LENGTH M#l 
'M °°. ORDER OF POLYNOMIAL . 
‘ROOTR RESULTANT VECTOR OF LENGTH M CONTAINING REAL ROOTS 

. ROOTI . RESULTANT VECTOR OF LENGTH M CONTAINING THE CORRE SPONDING 
“i. TMAGENARY ROOTS ; 3 
“TERY ERROR -CODE RETURNED WHERE : & 

: NO ERROR, ae a 


oPl= M LESS THAN LT 2 
. 2 = M GREATER THAN 36." 
. 3 = UNABLE TO DETERMINE ROOT. WITH 500 ITERATIONS ON 5 . 


“STARTING. VALUES 
HIGH ORDER COEFFICIENT 1S ZERO’ 


* 
1 


s. bad 


SUBROUTINE TAKEN FROM BM/SSP *. 
DIMENSION MOF ILI st F( 36 byRODTR(1),ROOTI( 1) 
IFYT=0 a 
N=M ; - 46 

TFER=0 \ 

LF ¢, XCOF(N#1) ) 20,25910. 

iy io IF'¢-N 1° 15,15,32 


PARAM MAMANMAANNRDAAANAADR|AAAGCAAA 


C ‘SET ERROR’ CODE TG 1 


15 I€R=l ; . OS 
20 RETURN ° a 


es : 
C SET ERROR CODE TO 4 : ‘ 
an i ; “ 


a0 


oo : 
ERIC 


25° TER=4 oo "i este Ss ? 
Go: 16.20 : : 2 eh oy 


c* ; fh “ s omer at is a 
7; SET ‘ERROR’ CODE TO 2 . or ee 2 a as 


"20. TER=2 pe eg 
60.59 20 . 2% » 3 
32. IF © N36) 95435530 . a ee a BS 
35  _NX=N : - a Se : 
NXX= Nel - & Sia e: Ky 
N2=1 . ; en 2 poe | LS a 
 KJIENtL eo Oe 


*  DO~GO L=19KJ1 eo Biz” = 
\ MT=KJteL +1 . 4 
40  COF(MT)=XCOF(L) is, © 2° eae 
C - - 
C SET INITIAL VALUES . ‘ 
c : : - 
45  XO0=.5001010E-2 | A ; ass 

YO=.1000101E~1 ote 


Cc ‘ 
‘C ‘ZERO INITIAL VALUE COUNTER 
Cc 


IN=0 SO . 2” & eee 


50 K=X0 ie 
u i 4 me 
GC. INCREMENT. oi ate AND COUNTER" ; : - 


XO=-10.4¥O0 2 


-  YO=-10e#X , ») ‘eS 4 an 
Cc ee . ae y ae 
C SET X ANDY TO CURRENT VALUE - = 
6 ° _ en es 
X=xXO_.. ‘ ns : ee ; — 
-.*¥2¥O * 4. °% ey es © | “3 
IN=IN¢1] | 
GO.TO 59 
55 IFIT=1 
XPR=X’ 
» ¥PR=Y 


: , . . 4 
© € EVALUATE POLYNOMIAL AND DERIVATIVES 


59_ ICT=0 

EV UX=0 6 r _ ee 5 
UY=0. . ; ‘ eA 
‘V=06 - P . 
YT=0. . . ¢ - 
: XT=l. * é 
U=COF(N+1) : 

. IF ( U } 654130765 

65 00 70 [=1yN a Z 

L=N-i+l _ ; 


BA Fut Text Provided by ERIC 
— 


ERIC 


C STEP ITERATION COUNTER — _ 


80 


85 


90 
C 


c. . SET ERROR cove TO 3 


eee 
95 


rc" 


105 


" ALO, 


115, 


120° 
aCe 


125 


440 


_ YT 2=X#YTHY*XT / 
— UXSUX4FI#XT#COFTL) 


_ YT=YT2 és 3 
SUMSQ=UX#UX+ UY#UY | i . 
5 . 


Go TO 20 


 ME=KSL-L+1 
FEMP=XCOF (MT) ; 


TTEMP=N - 


' 66,70 140 “ge 83 So . 


NEN-1 bom b* 


XT2= X#XTHYRYT 


USU+C OF (L)*XT2 | ‘ . ’ 
V=VHCOF(L)*YT2 _ i = 7 “os 
Fle! . 


UY=UY-FI*YT*COF(L) 
XT=XT2° . 


IF ( SUMSQ ). 75911097 

De CUE U TU RIL S UND 

X=X+DX 

DY =~ ( UFUY4 VEUX) /SUHSQ ; ‘ C ~ 
Y=Y+DyY : Pie , é 
IF { ABS ( DY) +ABS (DX) -1 .0E-5 ) 100, 80 80: a 


ICT=I1CT+#1 

IF ( ICT-500 ) 60,85,785 
IF ¢{ IFIT ) 100,790,100 
IF ¢ ane, ) AOU N72 


-TER=3 
DU 105 L=1yNXx 


XCOF (MT}=COF (L) “-S . eo ee : i 
COFULL)=TEMP a , & ; 


NENX . 
CNXSITEMP J tO 3 

"TFC LEIT 9120, $5 4120 ey 
LF ( LFIT ) 1155509115. . ee 
X=XPR OS 

Y=YPR. ae : 
TFIT=0° 7 ~ «a 
LHF OX ) 1229125,122 9, + 7 
IF ( ABSOY)~ABS (X)#1 .0E-4- ) 1354125, 125 
-ALPHA#X+X 

SUMS Q=X* X#\V# Y 

N=N-2) pe 


X=Q. £ . a * ‘ " i \ 
NX=NX-1 " 3 , 

NXX=NXX-1 
Y=0. ° 


SUMS Q=0. . se : 
AL PHA=X Ns : vies 


af 
— 


Ui=1 te ee . | 


* 


L2=2 - 
COF(L2)=COF(L2)-+ALPHASCOFELLD « 
Do.. 150 L= Zan 


ROOTI (NZ) =Y 
— -ROOTRIN2Z)=X 
‘+ NZ=N242 
IF -¢( SUMSQ ) 1605165160 
Y=e-¥ | 
SUMSO= Oe 
.GO TO 155 7 
IF ( N } 205-2024> 
END - 


1A Fut Toxt Provided by ERIC 


© SUBROUTINE SUBTST {IND ICH KPOP NDA SK 


Cc 
_f CALCULATION OF PARAMETER ESTIMATES FROM SUBTEST RESULTS 

C.. ; 

C PARAMETERS ESTIMATED: MEAN TEST SCORE : Z(1) 
‘{C ° , VARIANCE OF TEST SCORES - 202) 
ou" ~ °  3RD CENTRAL: MGMENT OF TEST ‘SCORES Z{3) - 
ee . 4TH CENTRAL MOMENT OF TEST SCORES. Z(4) 
VARIANCE COMPONENT FOR ITEMS ‘ z(5) 
ce © P (VARIANCE OF ITEM-OIFFICULTY INDICES) 

C, « és , VARIANCE COMPONENT FOR EXAMINEES Z6) 
-C VARIANCE COMPONENT FOR INTERACTION Z(7) 
_ oe % ESTIMATED RELIABILITY. OF SINGLE ITEM 2(8) 

c 


COMMON /BLOCK1/ SDATA( 102,11) NTEST % 
DEMENSION FMT(162) : 
*, DIMENSION XC 300}'s SPLUS(150) pSPLUS2 (150) sPIJ(11175 ) P6150) +208) 


ana 


NRDR=5 7 
IF ¢ NDISK, «67. O )- NROR=2 
C : 
Cl pate a 
~ WRITE? 1693) , " 
7 D0 500 L=LyNTEST . 
C eS 
READ 1597) NSPTsIPT 
. CALL ROFMT (FMT)NC)” . 
- -NCENC-1 © 
eS 
C IF THERE ARE MORE’ THAN 100 SUBTESTS, THE REMAINING SUBTEST RESULTS 
C ‘ARE IGNORED. ALSO, ANY SUBTEST HAVING MORE THAN 150 ITEMS IS 
C, IGNORED. 
aE 


KF (1 SLE. 100 AND. IPT «LE. 150 ) GO TO 75 
WRITE (695) I 
BE 50 J=LyNSPT 
50. READ (NRDR@FMT) (X(K) ¢K=1yIPT) 
" GU TO 500, 
78> CONTINUE - 


ary) 


Cc 7 eee 
: SDATA(1 99) =NSPT . 
SUATACI,1LOI=IPT : / 
WRITE (624) 1 WNSPT p IPT y (FMT (J) J= 1yNC) 


SY=0. | 

SYZ=0. 

SY 3=0. r ian 

SY4=0. = 
 SP=0-6 : { . ‘ 
aparUR 


‘ $P3=0. 
SP4=0. =. 
. IF (. NDICH-.EQ. 1 )-.60 To 105° 
NCOV=IPT*(IPT-211/2 
...D0 100 J=LsNCUV 
100 PIS(J)=0- 6 oo 
. 105 00.110 ‘s= =EyTPT . P 
"PP EI=0. 
'  - . SPLUS (J2=0. : 
110 SPLUS2(J)=0. 
DO 200 L=1L»NSPT_ - os ; ; 
READ (NRDR SFMT) “(X(K),K=1, IPT) 
- Y=0. cad « 
DO 120 K= 1,1PT 
T=X(K) 
“PC KD=PLK)4T 
- 120: Y=e¥eT 
NN=IPT-1 ° : 
BH 130 J=lyNN ‘ 
MM=J+1 : BAe 
DO 130 K=MM, IPT Z 
LSUB=(d-1) #1 PT#K- JH SE) /2 ; 
130) PES(LSUBJ=PIJ(LSUB)+X( 5) 4X(K)- 
»  ° PO 140 K=1,IPT 
SPLUS(K )=SPLUS(K)+Y¥#X¢K) /IPT 
740) SPLUSZ2(K)= SPLUS2(K)#¥#Y*X(K)/(IPT#IPT) . + tems ° 


SY =SY+Y _ ‘ee 
SY 2=SY24+¥*#2 i 


SY BES Y34+¥#43 
200 SY4=SY4+VH%4 


CALCULATION OF COMPONENTS. OF VARIANCE FOR SUBTEST MATRIX 


CXX=SY*SY/(NSPT*IPT) = ; t 
- SSE=S¥2/1 PT-CXX 
v=, 
-DG 300 J=1,IPT3 
BUD T=T+P (SFP CS) 
SS I=T/NSP T=C XX 
SS T=SY-CXX 
SSIE=SST-SSE-SS1_ 
SRR SSIE/( (LP 1-1, )*(NSPT~10)) 
2(5)=(SSI/(IPT-16)-ERRI/NSPT 
Z(5) #(SSE/(NSPT-159-ERR)/1PT 
24 T)=ERR 


FSTIMATION .OF RELIABILITY OF SINGLE ITEM 


FMSE=SSE/ (NS PT-1.3 

T=(NSPT~-1.6)* (IPT-1.) 
~FMP=(T-2.)/T 

Z2(8)=(FMSE-FMPKERR)/ (IPT*FMP¥ERR) 


i sec oe 7 ai 
eo« 
ot al 
¥ wpe aa i 
as beg *. * 
€ 2s 
: * 
a 


, - 98 350 J=1ly IPT 
PE3) =P (J) /NS PT ee 
_ TSPCury. : : 
“SP=SP4T 
| SP2=SP24+T#*2 
'  $P 3=SP3+T**3 = ; ; . . Cee 
350 SPAESE RTT SEG ; Bue - > Oe 


WRITE (6,6) (PEK) ) aK3h 2 LPT) . 


“ESTIMATE OF POPULATION MEAN AND VARIANCE WHEN ITEM “SCORED 
NON-D.ICHOTOMOUS LY . . 


one. 


IF (. NOICH «EQ. 0 ’ 60 TO 360 


oO 


2(.19=(KPOP /FLOAT( IPT) }#( SY/NSPT) 
2(2)=((SSE/( 1PT*(NSPT~1.)))- ((Le-FLOAT(IPT) /kPOP )#68R/ 181) 94 
1 KPOP*KPOP* (NS PT-1,.) /NSPT 

2(3)=0. . 

2(4)=0. a 


- GD TO 396 


' ese, Te gee : a 


BaAAGeEA 


CALCULATION OF G STATISTICS USED IN ESTIMATING MORENTS - 8 tg 


"ge 360 «G1 = =( SY/NSPT) #4 a 44 - a Sf C2 4 A 
-@. cant sv/nsr ty anes sY2/NSP1) eo Ce ee BA a, oe 
G3=SPHSPES P2 ee are ae ee oe ee 
G4 aCSY2/NS PT ) 2 Sk ee oe Seen a el ff 
G5=SP2*SP2 
G6=(SY/NSPT*(SY3/NSPT) ye & ee fae 
G7=SP¥SP3 | ; 6 8 ght MS oI a 8 2 ao 
. G8=SY4/NSPT _ . _— Ce. ea. at oe 
oe G9=SP4 - * ees ae. 
: G1 1=(SYZ/NSP T)*SP2 - fe ae 
G12=( SY/NSPT ) **3 => § 
G1 8=( SY/NSPT)*{SYZ/NSPT) 
GI9=SP*SP2.. 
G20=SY3/NSPT 
GZ1=SP3 
G23=(SY/NSPT )¥*200 © 
G24SSY2/NSPT a 3 Coon 
G25=SP2 . 4 2 ie eee 
G26=SY/NSPT 2 a 
G¥O=0. . 
Ce P 4 2 oe Gl 4=0. o ; ; ; 
ts G1 5506 6 2s 4 
ice 7 G16= Oe - ¢ ; ae i ss a os wat 
D0 375 J= L TPT - ~* 3 
G10=G104P (J) *SPLUS(J)/NSPT a Lee ob Bs 
G15=G15+P( J) *SPLUS2( J) /NSPT | ie ox yee ae 


“oO 375, G1 6=6 16+ P( J) #PL JS) #SPLUS (J) /NSPT 
in G22=G10*IPT 


, 36 7 os 


"ego ea eg 
. . . . Pa : ” . 


-G10=G10*SY/NSPT¥IPT aw x ue 
6155 GISF IP Te IPE. ae _ ¢ : ae 
“GLO=G1O#1PT | -,. @ 24 ——e 
OO OO NNSEPT=2 ie - & 
“!  ,.00.-390 ds Ly NN" _ ; ar 
ota) MM=J+1- 
. +: 00: 390° K= =HMy IPT. 
. LSUB=( U1) #1 PT4K~- 0 F( Sa /2. | 
390: CL4=G144P 14) ¥P CK) *PT SCLSUB) /NSPT ; ~ 
GL4=G1442645 P32 eS. _ oe, Ne 
'C, ESTIMATE OF POPULATION MEAN: AND VARI ANCE 
CI TEMS: SCORED. DIGHOTOMOUSLY } : 


ae 
es 
oe. : 
- 

7 


has exe 3 
B=A*( KPOP=1~)/(1PT-1.) 


* - CSBR(KPOP=2. UEP TH2 0 | 
D=C#LKPOPT 36 1/(18T~3. ) 


: ; “Lee 
“Ge 


~. “ZOLPEA*G26 © 
ADS we (626-6259 484 G24~626-G234625) 
spf aie =ZA2ZV#NS PT/(NSPTHLS ) 


ESTIMATE OF 3RD. AND’ 4TH GENTRAL MOMENTS © _ 
CLTEMS “SCORED DICHOTOMOUSLY) ar 7 


"2(3)=A*( 626-3 0*625 42. eGZLIE 


al B*(3.*G 24-3 .%626-3.*623+9. #G25+ 6.462246. ¥G19-G. *G21)+ 
2." C¥( 2046 26-3. Boe ct pcos a eeceeshGeet acter lac rteloit 
30 Gaye ln eetGse tes *Gi2). 

_ 


“24)= Ne626- = 44 #625 46'54G21-3.469)4 


‘43 B¥ (7 0*G26428 «*G25+7 « ¥624-4.%62 3-24. 4622-42 .#G21418- *G619+ 
2° 12 6¥G 16412. *G14421. #69-12.*G6 7-9 0*G5)+ ' 

a) C¥$122*G26-48 .*G25- 18. #62441 2.*G23+60» #62247 20462146046 20- 
4. 48.*G19-12 0¥G18-36. *Gi6—-12 «#615 -24 6 G1 446 o¥G1dt60 *G1Li +240 G10 
§ ~360*G9 £364 4G7418 o*G65-1 80463) + 
6 D* ( -6 6* 626.411 0624-66 #620468 +24 0*G25-8.*623- = 36.%G22412.%6 1a4 
qT 126 ¥G15~454G66-360%G214300%G1 94246 #616412 6¥*614-66%*G612-66*G11-_ 
8 24 2610 460%62418 - ¥69-24.6*G67- 9 *G5418.*G 3-3. *G61) 

396 00° 400 J=1ly,8 ; ; . 


400 SDATA(I,JS)=Z2(05) 


500 CONTINUE . 
-" IF ( NIEST .GT. 100, } NTEST=100 


WRITE (6,2) a : 
00 550 I=1yNTEST ; ; * 


cy . 7 + 
, Zz ” 


FORMAT (/7725X521HESTIMATES ‘OF HOMENTS 133%, ge ge OF ie 
‘L22HCOMPONENTS OF VARIANCE/8X,__ Oe : . 


2B Tap mmm nny enn teen bpm BK 
J ByS Hema nn nnn an Son ne on beens st 1 


- 48H SUBTEST s5Xs4HMEANs 9X pBHVAR TANCE ¢ '10Xy 2HU3 91 3X 2HUG 
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